Spin-triplet pairing in large nuclei 
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The nuclear pairing condensate is expected to change character from spin-singlet to spin-triplet 
when the nucleus is very large and the neutron and proton numbers Z, N are equal. We investigate 
the transition between these two phases within the framework of the Hartree-Fock-Bogoliubov 
equations, using a zero-range interaction to generate the pairing. We confirm that extremely large 
nucleus would indeed favor triplet pairing condensates, with the Hamiltonian parameters taken 
from known systematics. The favored phase is found to depend on the specific orbitals at the Fermi 
energy. The smallest nuclei with a well-developed spin-triplet condensate are in the mass region 
^ ~ 130 - 140. 



INTRODUCTION 



The neutron-proton interaction in the spin-one chan- 
nel is attractive and stronger than the identical-nucleon 
interaction. Nevertheless, the nuclei observed in nature 
favor identical-particle pairing. A trivial answer to this 
puzzle is that most nuclei have different numbers of pro- 
tons and neutrons, and the isospin polarization discour- 
ages pairing in the isospin-zero (T=0) channel. But even 
in nuclei with equal numbers of neutron and protons the 
J = 1,T = neutron-proton pairing does not compete 
well with ordinary J = 0,T — 1 pairing, as evidenced by 
the ground-state spins and binding systematics of odd- 
odd nuclei P,[3. This puzzle presents a serious challenge 
to nuclear theory that must be resolved if the theory is 
to be used confidently in contexts where experimental 
information is not available, such as astrophysical envi- 
ronments. 



It has recently been suggested that the explanation 
of the suppressed spin-triplet pairing^s the presence of 
the strong nuclear spin-orbit field [IHS, H, 11] ■ The spin- 
orbit field interferes with the pairing in both channels, 
but it suppresses it more strongly in the J — 1 neutron- 
proton channel than in the J = identical-particle chan- 
nel. Since the spin-orbit splitting is a surface effect, one 
might expect that the pairing would change character to 
the J — I condensate in extremely large nuclei. In the 
limit of large nuclear size, the surface-to-volume ratio 
would be small and the spin-orbit field would be ineffec- 
tive at controlling the single-particle spectrum. In this 
paper we apply the Hartree-Fock-Bogoliubov (HFB) the- 
ory with a simplified Hamiltonian to examine the proper- 
ties of the pairing condensate as one approach this limit. 
We shall find that the transition to J = 1 pairing does 
occur in large nuclei, but is only strong for nuclei that 
are somewhat beyond the limit of stability with respect 
to proton emission. 



THE HAMILTONIAN 

To keep the calculations as simple as possible, we con- 
sider a Hamiltonian consisting of a phcnomenological 
single-particle Hamiltonian together with a two-particle 
interaction, 

H = ^{i\Hsp\j)alaj + ^ (ij>|/c/)aja]a,afe (1) 

i i>j,k>l 

Here the single-particle Hamiltonian Hsp is taken with 
the well-known Woods-Saxon form. 



P 

2m 



+ Vwsfir) + I ■ sV,o 



1 d!{r) 
r dr 



(2) 



with f{r) = 1/(1-1- exp((r — R)/a) and R the radius 
of the potential. We use parameter values Vws = —50 
MeV, a = 0.67 fm which are close to standard values in 
the literature 0, Eq. 2.181]. We do not expect that the 
competition between the different kinds of pairing will 
depend sensitively on Vws and a. 

The spin-orbit strength parameter Vso will be impor- 
tant in the competition, and its determination requires 
some discussion. Unfortunately, first-principles studies of 
the nuclear Hamiltonian are not at a point where quan- 
titatively reliable spin-orbit fields can be calculated 0]. 
However, given the Woods-Saxon form in Eq. ([2|), one 
can constrain the parameter quite well from spectroscopic 
data. A good source is the energies of single-particle exci- 
tations in nuclei adjacent to closed shells. For the present 
study we fit the experimental spin-orbit splitting of the 
f-!/2 and /5/2 in the nucleus''^^ Ca The measured value 
is 7.0 MeV, and the fitted spin-orbit strength is Vso = 33 
MeV-fm^ . We estimate an uncertainty of not more than 
10% in this value, based on comparison with other ways 
to estimate it. 

We only treat nuclei with equal numbers of neutrons 
and protons in this work, enforcing the condition by tak- 
ing equal chemical potential in the HFB equations. The 
value of the chemical potential is fixed by demanding 
that the density goes to the standard value p = 0.16 
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frn"^ in the nuclear-matter limit. With the chosen single- 
particle Hamiltonian, the corresponding chemical poten- 
tial is A = —13 McV. Note that our Hamiltonian neglects 
the Coulomb interaction. This allows us to approach the 
nuclear matter limit but the calculated nuclei extend to 
the unphysical region beyond the proton dripline. 



Interaction 

The two-particle matrix elements are determined by a 
contact interaction of the general form 



given shell-model spaces, it is still necessary to choose a 
suitable space cutoff parameter Eb- An argument can be 
made identifying Eb with the harmonic oscillator shell 
energy fiw, which would put it in the range of 11 — 14 
MeV.for the strengths in Table I. For our study here, we 
have somewhat arbitrarily taken the strength of the spin- 
singlet and spin-triplet interactions to be Vg — 300 MeV- 
fm'^ and Vt = 450 MeV-fm'^, respectively, with a cutoff 
of Eb = 10 MeV. The ratio is 1.5, giving the spin-triplet 
somewhat less of an enhancement than the fits suggests. 
The calculations will then be on the conservative side in 
testing for a possible phase transition. 



{tj\v\kl) = ^z;„(u|<5(3n^"-^"')^L=o^'a|fcO- (3) 



Here Pl=o is a projection operator on total orbital an- 
gular momentum zero. The Pa is a projection op- 
erator on the six spin-isospin channels {Ss^,Tt^) — 
(Oo, li), (Oo, lo), (Oo, l-i), (lo, Oo), (li, Oo), (l-i", Oo)' The 
interaction strengths Va are of course independent of the 
^-projections of spin and isospin, and we shall call the 
two independent strengths and Vt for spin-singlet and 
spin-triplet, respectively. However, we will need to keep 
explicit the six amplitudes with different z-projections 
when constructing the HFB condensate. It should also 
be noted that the contact interaction has to be regulated 
in some way to produce a finite condensate. This will 
be implemented by including only orbitals whose orbital 
single-particle energies ei are within a range Eb/'2 of the 
Fermi energy, X - Eb/2 < ei < X + Eb/2. 

Our choices for the parameters will be guided by the 
relevant shell model matrix elements of phenomenolog- 
ical shell-model Hamiltonians. In this way we hope to 
avoid some of the ambiguity associated with the form of 
the pairing interaction (see concluding section). Specif- 
ically, we take the USDB Hamiltonian fitted to sd-shell 
nuclei llOl and the GXlA Hamiltonian fitted to fp shell 
nuclei [llj. We consider all the shell- model matrix el- 
ements with total spin and isospin couplings (J, T) = 
(0, 1), (1, 0), and make a least-squares fits to each set us- 
ing the contact interaction and harmonic oscillator or- 
bitals. In the fit, we give equal weight to all matrix el- 
ements in the m-scheme representation. The results for 
the ratio of interactions strengths vt/vs, shown in Table 
I, are in the range 1.6-1.7. Another fitted Hamiltonian 
to /p-shell spectra [Ij gives a very similar ratio (1.56). 
To get the absolute values of the strengths, we also need 
to specify the oscillator parameter controlling the size of 
the orbitals. This is taken as fiw = 14 MeV and 11 MeV 
for the sd-shell fit and the /p-shell fit respectively. The 
deduced strength parameters Vg^Vt are also given in the 
table. One sees that the interactions derived from USDB 
and GXlA are in good accord with each other. Unac- 
countably, the interaction from ref. 0] has much lower 
strengths. For application to the HFB theory outside 



TABLE I: Strengths of triplet and singlet interactions from 
shell-model fits and their ratios. See text for details. 



Source 


Va Vt 

MeV-fm^' MeV-fm^ 


Ratio 


sd-shell \10} 
/p-shell [nl 


280. 465. 
291. 475. 


1.65 
1.63 



HFB 

The physical quantity to be calculated is the energy of 
the system in a state constructed by making a Bogoliubov 
transformation on the noninteracting ground state. The 
expression for the energy is 

E ^^Pij < i\Hsp\j > + ^Y1 

ij i>j i'>j' 

(4) 

where pij and Kij are the ordinary and anomalous densi- 
ties associated with the transformation. There is also a 
term in the interaction energy which we ignore along 
with all other Hartree-Fock contributions. Another con- 
venient quantity for studying the pairing phases is the 
condensation energy, also called the correlation energy 



Er, 



Eo-E. 



(5) 



where Eq is the energy of the ground state in the absence 
of a pairing condensate. The HFB ground state will max- 
imize the correlation energy. By the Bloch-Messiah the- 
orem 12|, p. 611-612], the solution of the HFB equations 
could also be found by the BCS minimization with the 
variational wave function \BCS) = Yiai^a + VaaaCisl), 
provided one knows the orbitals a and the pair corre- 
spondence a ^ a. For ordinary spin-zero pairing in nu- 
clei with even numbers of nucleons, we know that the 
pairs are the partner orbitals under the time-reversal op- 
erator. If the orbital basis does not have more than one 
orbital of given orbital angular momentum i, the partners 
are related by the time-reversal operator, \a) = \£sj —jz) 
for I a) = \isjjz). Thus, for narrow energy-band windows 
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Eb and not too large nuclei, the BCS theory is adequate. 
This will provide one check on the more complicated HFB 
calculations. 

For spin-triplet pairing in the presence of a spin-orbit 
field, the HFB treatment is unavoidable because the op- 
timum orbital basis \a) depends on the relative strengths 
of the pairing and spin-orbit fields. However, there is one 
case where the BCS treatment is valid. That is for a sin- 
gle j-shell and pairing in the 8^ = channel. In that case 
the pairing field demands that jzict) — ~jz{a). This is 
uniquely satisfied by the same partnering as in the S = 
pairing channel. We shall also use this case for testing 
the full HFB code. 

Another important physical quantity is the quasipar- 
ticle energy Eq, defined as the eigenvalue of the HFB 
matrix equation, Eq. ([6]). It has the interpretation as 
the removal or addition energy for a particle in system 
with an even number of nucleons. 



To efficiently calculate the pairing field elements Ay in 
the HFB matrix, we save and store the anomalous den- 
sities Ka(r) on a radial grid. The densities Pij and 
are calculated in the orbital representation IJ, p. 251, 
eq. (7.23)] from the selected columns of the transfor- 
mation matrix that diagonalizes the HFB Hamiltonian. 
Under the assumption of spherical spatial symmetry, the 
HFB matrix is block-diagonal with respect to the orbital 
quantum number i; the different ^-blocks are calculated 
separately in our codes. If one restricts the anomalous 
density to channels having Sz = 0, the matrix can fur- 
ther be decomposed into blocks of fixed \jz \ = + Sz|- 
Further details are given in Appendix A. 



RESULTS 
An example: ^*Cr 



The 5'I7(4) limit 

Before presenting calculations with physical values of 
the Hamiltonian parameters, we examine the theory in 
the limit of SU (4) symmetry. The symmetry requires 
that the two pairing strengths be equal and that the 
potential be purely central. The HFB ground state of 
this Hamiltonian is highly degenerate, with many dis- 
tinct ways to form the pairs. For example, one solution 
of the HFB equations is the state having both NN and 
PP pairing, with independent condensates for both. But 
we can equally well form a condensate pairing up-spin 
neutrons with up-spin protons and another condensate 
pairing the down-spin nucleons with each other. Both 
states have the same energy in the SU{4:) limit. This 
will be shown explicitly the example presented in Sect. 
IV. More generally, there is a continuous group of trans- 
formations that leaves the condensation energy invariant. 
The degeneracy of the HFB ground state can only be bro- 
ken by treatments of the wave function beyond the HFB 
approximation, e.g. [23| . 

CALCULATIONAL DETAILS 

We set up HFB equations in a basis of states |i) con- 
structed from the eigenvectors of the central Hamiltonian 

/2m + VwsF{r). We assume spherical symmetry and 
represent the orbital wave functions by their quantum 
numbers i = {n,£,£z, Sz) and radial wave functions 4>i{r) 
on a uniform radial mesh. The basis is truncated by keep- 
ing only states whose single-particle energies ££ (without 
spin-orbit) are within a range Eb /2 of the Fermi energy, 
|££ — cfI < EbI^. The spin-orbit interaction is treated 
by including its orbital interaction matrix elements in the 
HFB Hamiltonian. 



Before applying the theory to very large nuclei, we 
investigated its performance in an experimentally well- 
studied region, namely the mass region corresponding 
to an open /7/2 shell. The middle nucleus filling the 
/7/2 shell is "^^Cr, with 4 neutrons and 4 protons outside 
the completely filled shells of Is, Ip, 2s, and \d orbitals. 
To model this nucleus, we take the radius parameter 
R = 4.62 fm corresponding to the phenomenological The 
single-particle Hamiltonian gives a spectrum that corre- 
sponds very well with the shell assignments with the /7 /2 
shell at the Fermi energy. The HFB equations are solved 
adjusting the chemical potential to get the correct parti- 
cle number, A = 48. Some results for the HFB theory are 
presented in Table II, showing how the gaps and correla- 
tion energies depend on the Hamiltonian properties. We 
first demonstrate the degeneracy of HFB solutions when 
the Hamiltonian has S'i7(4) symmetry in the spin and 
isospin degrees of freedom. The required Hamiltonian 
has spin-field field set to zero and equal pairing strengths 
Vs and vt- The calculated properties of the condensate 
are shown in the first four rows of Table II. One sees that 
the calculated quasiparicle energy is independent of the 
quantum numbers of the condensate and is the same for 
all orbitals. The correlation energy is also independent of 
the choice of condensate, demonstrating the degeneracy 
of the S'C/(4) HFB ground state. 

Going now to a more realistic Hamiltonian, rows 5 and 
6 in the Table show the maximum effect of the spin-orbit 
field, taking only the /7/2 shell for the single-particle 
space. Row 5 shows the results for an ordinary spin- 
singlet condensate, pairing neutrons with neutrons and 
protons with protons. The quasiparticle energy Eq is the 
same for all 8 orbitals in the shell. The value, Eq \ 
MeV, is in fact close to the average experimental odd- 
even mass difference. Row 7 shows the results for spin- 
triplet pairing. The S = \ condensates define a direction 
in space and the gaps are no longer independent of the jz 
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quantum number of the orbitals. These gaps are plotted 
as a function of jz in Fig. 1. One sees that the gap ap- 
proaches zero as jz becomes small. This behavior would 
be called gapless superfluidity in a large system. The 
correlation energies of spin-singlet and spin-triplet pair- 
ing can be compared in the last column of the Table. One 
sees that it smaller for spin-triplet than for spin-singlet. 
Thus, the ground state should exhibit ordinary pairing, 
as expected. 

We next augment the single-particle space by adding 
the shell, including the physical spin-orbit splitting. 
The results for four choices of condensate are presented. 
First of all, one sees that the correlation energies are 
larger than in the pure /7/2 Hamiltonian. They could 
not be smaller because the HFB theory is variational: cn- 
largening the space can only lower the energy. The first 
two rows here show the results for spin-singlet pairing, 
differing on the isospin coupling. The energy are exactly 
the same due to the isospin invariance of the Hamilto- 
nian. A similar degeneracy is found for the spin-triplet 
pairing; here the condensate energy does not depend on 
Sz- We shall exploit this invariance later by limiting our 
trial condensates to be either (0, Iq) or (lo,0). Then jz 
is conserved and the HFB matrix can be diagonalized in 
small blocks (See Appendix A). The final entry in the 
Table is for the Hamiltonian in the full /p-shell. This 
is in fact the truncation that result from the Eb = 10 
MeV cutoff around the Fermi energy. We see that the 
correlation energies are larger, as expected. The spin- 
singlet pairing is still the stronger one, so the Hamilto- 
nian with the Eb passes the test of agreeing with known 
phenomenology. 

Systematics 

In this section, we compare correlation energies for the 
triplet and singlet pairing channels as function of nuclear 
size. We showed in the previous section that the corre- 
lation energy does not depend on the z— components of 
the spin or isospin of the condensate, so we may limit our 
exploration of possible condensates to the Sz = ^ chan- 
nels. As mentioned above, the HFB matrix reduces to 
blocks of fixed \ jz\ for these channels. We find the self- 
consistent solutions of the HFB equations iterating from 
a starting point in which there is a finite condensate in 
some channel. 

For the range of nuclear sizes that we consider, the 
pairing gap is comparable to or smaller than the energy 
spacing of the shell orbitals. Under these conditions the 
strength of the pairing condensate will be quite sensitive 
to the Fermi level. Since we are interested here in nuclei 
that have well- formed pairing condensates, we shall only 
calculate systems where the chemical potential (A = —13 
MeV) coincides with the orbital energy of some j-shell. 
This gives about 50 cases for nuclei with radii in the 



TABLE II: Pairing gaps and correlation energies for Cr in 
the HFB theory. D is the dimension of the A or k matrix; the 
dimension of the HFB matrix is twice that. The spin-singlet 
interaction strength is — 300 MeV-fm"^. The spin-triplet 
interaction strength is vt = 450 MeV-fm'' except for the rows 
with the 517(4) results; there Vt = 300 MeV-fm^. The con- 
densates are labeled by the spin-isospin quantum numbers, 
dropping the z-quantum number when S = or T = 0. The 
column marked Eq gives the value or the range of values of 
the quasiparticle energy of orbitals at the Fermi energy. The 
correlation energy Ecorr in the last column is defined in eq. 



Space 


D 


condensate 


Eq 


Ecorr 






(Ss^.TtJ 


(MeV) 


(MeV) 


/-shell 


28 


(0,li),(0,l_i) 


1.80 


10.27 


SU{4) 




(0, lo) 


1.80 


10.27 






(li,0), (1-1,0) 


1.80 


10.27 






(lo,0) 


1.80 


10.27 


/7/2 


16 


(0,li),(0,l-i) 


1.03 


4.11 






(0, lo) 


1.03 


4.11 






(lo,0) 


0.12-0.87 


2.01 


/7/21 


28 


(0,li),(0,l_i) 


1.16 


4.64 






(0, lo) 


1.16 


4.64 






(li,0), (1-1,0) 


0.20-1.34 


2.98 






(lo,0) 


0.20-1.34 


2.98 


full fp 


28 + 12 


(0,li),(0,l_i) 


1.23 


4.92 






(0, lo) 


1.23 


4.92 






(li,0),(l_i,0) 


0.26-1.58 


3.45 






(lo,0) 


0.26-1.58 


3.45 



range of 4.0-12.0 fm. These correspond to mass numbers 
in the range A = 2h — 1000. A table of the results for 
the two correlation energies and their ratio is given in 
Appendix B. For the lighter nuclei, correlation energies 
are of the order of several MeV, except for several nuclei 
with j = 1/2 shells at the Fermi level, for which the spin- 
singlet correlation energy can be less than 1 MeV. In the 
heavy nuclei, the correlation energy can be several tens 
of MeV. 

Fig. 2 shows a plot of the ratio of spin-triplet to spin- 
singlet correlation energies as a function of the mass A. 
The results show that there is a trend to favor spin-triplet 
pairing for large nuclei, as was argued in the Introduc- 
tion. However, the ratio is by no means monotonic as a 
function of R. Spin-singlet is favored for R < 6.5 fm and 
spin-triplet for R > 9 fm, but nuclei in between could 
have either ground state. The lightest nuclei predicted 
to have spin-triplet condensates are indicated by their 
mass and element designation. They are ^°P, ^^Sr, and 
^^^Er. The two lighter ones are in the physical region, 
and the spectroscopic properties of one of them, '^"P, are 
well-known. The calculated correlation energy for this 
nucleus is 1.4 MeV, which is not large enough to be con- 
sidered a strong condensate. Nevertheless, the observed 
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FIG. 1: Quasiparticle energies in **Cr for the /-shell space. 
Red circles: spin-singlet; blue squares: spin-triplet with con- 
densate in the Sz = channel. Lines are drawn to guide the 
eye. 



ground state spin and parity (J'^ = 1+) agrees with the 
quantum numbers of the spin-triplet condensate. The 
middle nucleus, ^^Sr has a somewhat smaller condensa- 
tion energy, Ecor ~ 1 MeV, so any effects of the pairing 
would be quite weaker. The nuclei in this mass region are 
accessible with current-generation accelerators and their 
experimental spectroscopic properties are under active 
study. The third nucleus, ^"^^Er, just beyond the physical 
region, is the smallest nucleus having both strong pair- 
ing {Ec ~ 13 MeV) and a spin-triplet condensate. There 
are two nearly degenerate shells at the Fermi energy, the 
2si/2 and lci3/2 shells. They both have radial nodes and 
weak spin-orbit splitting from their j-shell partners. 

A table of the correlation energies is given in Appendix 
B. The table also shows the angular momentum quantum 
numbers at the Fermi energy for each case calculated. 



CONCLUSION 

This study offers a resolution of the conundrum, why 
is it that ordinary nuclei only exhibit spin-singlet pairing 
when nuclear matter calculations find that spin-triplet 
pairing is stronger? The answer, that ordinary nuclei are 
too small to be out of the influence of the surface spin- 
orbit field, is seen by calculating the pairing energies over 
a large range of sizes. For our Hamiltonian, the spin- 
triplet pairing dominates for all calculated nuclei whose 
radii are larger than ~ 9 fm, and for a number of lighter 
nuclei having low-Z shells at the Fermi energy. The can- 
didate for the smallest nuclei with a well-developed con- 
densate, A ^ 130— 140, is tantalizingly close to region of 
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FIG. 2: Ratio of spin-triplet to spin-singlet correlation ener- 
gies as a function of mass number A. Nuclei with spin-singlet 
and spin-triplet condensates are shown as red circles and blue 
squares, respectively. The vertical line at A ~ 120 shows the 
dividing line between nuclei that are bound (left) and nuclei 
that unstable with respect to proton emission, according to 
the mass table of ref. 1 1311. 



physical nuclei defined by the single-nucleon driplines. 

However, due to the simplifications made to the Hamil- 
tonian and to other uncertainties we cannot be quan- 
titative about the transition point. In particular, im- 
portant physics may be lost in ignoring the momentum- 
and density-dependence of the pairing interaction. In the 
theory of nuclear matter the spin-singlet pairing derived 
from Brueckner theory is very weak at nuclear matter 
density [li. See Fig. 8][l3; see also the contradicting 
claim in ref. [l6| . The pair gap becomes larger at subnu- 
clear densities, suggesting that surface effects are crucial 
for that channel. In many studies of nuclear structure, 
the authors assume that there is some surface enhance- 
ment of the pairing, see e.g. 



17, 18, h 



Beyond the 

effects that are treated in the Brueckner theory, there is 
likely to be a significant induced interaction associated 
with collective surface vibrations [l^l . If we included ei- 
ther mechanisms enhancing the surface region, the pair- 
ing energy would decrease as the nucleus becomes larger. 
This of course could effect the crossover from singlet to 
triplet condensates. Also one needs to take into account 
renormalization effects associated with the truncation of 
the orbital space. One finds that the renormalized inter- 
action is actually suppressed in the surface region 21 1. 
Clearly, there is much left to do to improve the interac- 
tion. 

It would also be interesting to explore the effects of 
nuclear deformation. Breaking the spherical symmetry 
would effectively reduce the single-particle level density 
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at mid-shell filling, weakening the pairing condensates. 
How this affects the competition between the two kinds 
of pairing is not known. 

There are a number of other questions that could be 
addressed within the context of the simplified Hamilto- 
nian that we have employed. Can two kinds of conden- 
sates coexist in the same system? At first sight this seems 
unlikely because the bulk pair correlation length is large 
compare to the sizes of the nuclei we consider. How- 
ever, it has been claimed that the Cooper pairs have 
a small size in the nuclear surface [22|. This might al- 
low the existence of a mixed phase with a strong spin- 
triplet condensate in the interior and a spin-singlet con- 
densate on the surface. Another interesting question is 
the fragility of spin-triplet pairing is when there is a neu- 
tron excess. For our midmass triplet-pairing candidate 
nucleus, = Z ~ 65, the condensate seems quite robust 
with a calculated correlation energy 6 MeV larger than 
the singlet energy. It might be that the system could sup- 
port a few excess neutrons and still preserve the triplet 
condensate. 

We have not addressed the question of how the triplet 
pairing condensate could be observed, beyond noting that 
the character of the condensate controls quantum num- 
bers of the ground state in an odd-odd system. In fact 
several J = 1,T = ground states are sprinkled among 
the 12 odd-odd N — Z nuclei known experimentally. The 
/3-decay spectrum may also be affected by the pairing (23j. 
Finally, two-particle transfer reactions would be affected 
by the coherence of the pairing fields. 

The authors thank M. Horoi for supplying us with the 
GXlA interaction. We also thank A. Bulgac, T. Duguet, 
B. Spivak and S. Baroni for discussions, and J. Vinson 
for help in an early stage of this work. This work was 
supported by DOE Grant DE-FG02-00ER41132. 



Appendix A 

The Hartree-Fock-Bogoliubov equations have the gen- 
eral structure 



Al 



-A* 



-Hf, 



Al 



(6) 



Here Hsp and A are matrices in some basis of orbital wave 
functions and Ua,Va are column vectors. The ordinary 
and anomalous density matrix elements are given by 



Pi] 



(7) 



the nucleus is spherical with spherically symmetric fields 
except for the spin variables, using as a basis of the or- 
bital wave functions the eigenfunctions of the central po- 
tential which may be labeled by the quantum numbers 

To efficiently calculate the pairing field elements Ay 
in the HFB matrix, we save and store the anomalous 
densities Ka{i') on a radial grid. The densities are defined 

'^"(^) " ^ XI (t>i{r)4>j{r){a\ii)K,j (8) 



where 



{ij\a) = V2((l/2s,a/2s2j|5,5,„) X (9) 



(1/2 ^.(z) l/2UXj)\T{a)TAa)) {-t~'-&u,iM.a.y 

and the radial wave functions 0i are normalized as 
r^dr|0i(r)p = 1. The elements of the HFB A matrix 
are calculated as 

/>oo 

Aij ='^{ij\a)va r'^dr(j3i{r)(j)j{r)Ka{r). (10) 

where Va is the strength of the contact interaction defined 
in eq. (3). 

Due to the symmetry, the condensation energies can 
be calculated under the further assumption that the to- 
tal azimuthal angular momentum \ jz\ — \elz-'r Sz\ is con- 
served. The HFB matrix then decomposes into 2 x 2 or 
4x4 dimensional blocks with respect to spin and orbital 
angular momentum and t — z. A typical subblock for 
a fixed isospin of the pairing field has the form 



hi,j^ - Al 



A. 



hf 





3^ 





Al 





-A, 




Al 



A, 






-hf,-j^ + Al. 
(11) 

Assuming that there is only one radial state in the Eb 
window, the h and A matrices have dimension 1 or 2 
depending on whether \jz \ = ^ + 1/2 or not. In the latter 
case, the hgj^ matrix is given by 



Si + jjz - l/2)w i/2 
V(^+ 1/2)2 -j>,/2 



V(^+ 1/2)2 -j>,/2- 
£e - Uz + l/2)mwe/2 
(12) 

where Si is the orbital energy calculated without the spin- 
orbit field and Wi is the radial matrix element of the 
spin-orbit field. The matrix A^ is given by 



where the sums are restricted to eigenvectors satisfying 
Ea > 0. We construct the HFB matrix assuming that 



di 
±dt 



(13) 



Here the upper and lower signs apply to the (ST) = (01) 
and (10) channels, respectively. The pairing field de is 



identical to Aij in eq. (fTU]) with eg. 



\n£izSztz) and 



7 



|j) = \n£ —£z —Sz —tz). Note that d( is the same for ah 
the HFB subblocks of given i. The 2x2 HBF matrix in 
the = £ + 1/2 case can be obtained by ehminating 
the orbitals with unphysical m values from the basis. 

To solve the HFB equations, we start with a trial field 
for the anomalous density Ka{r), typically a constant 
value independent of r in some channel a and zero in 
the other channels. The solution of the HFB equation 
with the trial field is used to calculate a new field, and 
the process is iterated to convergence. 

A computer code to solve the HFB equations based on 
these formulas is available by E-mail from the authors, 
bertschOu . Washington . edu. 

Appendix B 

Columns in the table refer to: 
R, the radius parameter of the Woods-Saxon potential 
eq. (2) 

A, the nuclear mass number 

Ec{s), the spin-singlet condensation energy in MeV 
Ec{t), the spin-triplet condensation energy 
Rts — Ec{t)/Ec{s), their ratio 

£,j, orbital and total angular momentum of the single- 
particle orbital at the Fermi energy. 
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